Thermodynamics and Spin Tunneling Dynamics in Ferric Wheels with Excess Spin 
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We study theoretically the thermodynamic properties and spin dynamics of a class of magnetic 
rings closely related to ferric wheels, antiferromagnetic ring systems, in which one of the Fe (III) 
ions has been replaced by a dopant ion to create an excess spin. Using a coherent-state spin path 
integral formalism, we derive an effective action for the system in the presence of a magnetic field. 
We calculate the functional dependence of the magnetization and tunnel splitting on the magnetic 
field and show that the parameters of the spin Hamiltonian can be inferred from the magnetization 
curve. We study the spin dynamics in these systems and show that quantum tunneling of the Neel 
vector also results in tunneling of the total magnetization. Hence, the spin correlation function 
shows a signature of Neel vector tunneling, and electron spin resonance (ESR) techniques or AC 
susceptibility measurements can be used to measure both the tunneling and the decoherence rate. 
We compare our results with exact diagonalization studies on small ring systems. Our results can 
be easily generalized to a wide class of nanomagnets, such as ferritin. 
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I. INTRODUCTION 

Nanomagnets and molecular clusters are systems in 
which macroscopic quantum phenomena may be ob- 
served irt-the form of quantum tunneling of the magneti- 
zationjjij Two scenarios must be carefully distinguished: 
incoherent macroscopic quantum tunneling (MQT) and 
macroscopic quantum coherence (MQC). In the latter 
case, tunneling between energetically degenerate spin 
configurations takes place at a rate A/h large compared 
to the spin decoherence rate F. In ferromagnetic molecu- 
lar clusters such as Feg and Mni2 the aroimd state tunnel 
splitting A is small compared to ?iFfl Q However, 
significantly larger in antiferromagnetic (AF) systemd 
which are promising candidates for the observation of 
MQC in the form of coherent tunneling of the Neel vec- 
tor n. 

ThgJerric wheels (FWs) LiiFeg, Na:Fe6, CsiFcg, and 
Feioullj are a particularly interesting class of molecu- 
lar magnets. The s = 5/2 Fe (III) ions are arranged on 
a ring, with an AF nearest-neighbor exchange coupling 
J > 0, and a weak, easy-axis anisotropy {kz) directed 
along the ring axis e^. For h^ — 0, the classical ground- 
state spin configuration has alternating (Neel) order with 
the spins pointing along ie^. The two states with the 
Neel vector n along ±e^ [Fig. |l|, labeled | t) and | |), 
are energetically degenerate and separated by an energy 
barrier of height Nk^s^ . However, | t) and | |) are not 
energy eigenstates. Rather, the low-energy sector of the 
FW consists of a ground state \g) = (| T) + I i))/'s/2 and 
a first excited state |e) = (| T) — I i))/V2, separated in 
energy by A. For weak tunneling, \g) and |e) are ener- 
getically well separated from all other energy eigenstates. 

For rEaw , the tunnel splitting A can be as large as 
2.18KO'tj An estimate for the electron spin decoher- 
ence rate F in FWs can be obtained from the typical 



energy scales of the various interactions of the electron 
spins. These include nuclear dipolar interactions with 
^H nuclei (O.lmK [Ref. [16[), hyperfine interactions with 
^^Fe (ImK [Ref. |l^), and interring electron spin dipolar 
interactions (of order 10-50mK), and possible interring 
superexchange processes, which are difficult to estimate 
but may be of order lOOniK. However, the tunnel split- 
ting is sufficiently large that FWs are most promising 
candidates for coherent tunneling of n even in the pres- 
ence of intrinsic or extrinsic sources of decoherence on 
energy scales up to 0.5K. 




FIG. 1. The two degenerate classical ground-state spin con- 
figurations of the ferric wheel Feio. 

Although quantum tunneling effects in antiferminag- 
nets are more pronounced than in ferromagnetsjj'u the 
detection of quantum behavior is experimentally more 
challenging. The reason for this is that magnetization 
and susceptibility measurements probe only the total 
spin of the molecule which, by symmetry, remains un- 
altered upon tunneling of the Neel vector. This problem 
is resolved in ferrimagncts and antiferromagnets with un- 
compensated sublattice spins, in which A is still large and 



tunneling of n leads to lar, 
current (AC) susceptibility,!! 




nals in the alternating- 
ivided that magnetic 



cases J' / J » 1 and J' / J <C 1. The relevance o f the 



fields are smallHiTc^ Recent worked indicates that ferri- 



magnets exhibit a wealth of interesting tunnel scenarios 
also in finite magnetic fields. 

In this paper we study magnetic rings closely related 
to the ferric wheels, in which one of the Fe (III) ions 
has been replaced by a Ga or Cr ion with spin s' ~ 
or s' = 3/2, respectively, to create an excess spin. Such 
systems have been synthesized recently.EEl Starting from 
a microscopic model Hamiltonian, 

N-l 
H = J^Si- Sj+i + J'(si • §2 + Si • Sat) 
j=2 

N N 

+ Sh.^S,-(fc;s?,, + fc.^s2j, (1) 

1=1 1=2 

which also accounts for modified exchange (J') and 
anisotropy (fc^) constants due to doping, we calculate 
various thermodynamic quantities and spin correlation 
functions. In Eq. ([|), A^ = 10, 8, or 6, h = qhb^/K with 
B the external magnetic field, and g = 2 \s the electron 
spin g-factor. As we will show, the excess spin 6s = s' — s 
is strongly coupled to n and hence is expected to modify 
both the thermodynamic properties and the spin dynam- 
ics of the FW. In contrast—for an impurity spin coupled 
weakly to the Neel vector,c3 the thermodynamic proper- 
ties of the wheel remain essentially unaltered. For the 
modified FWs discussed in the present work, Neel vector 
tunneling also leads to oscillations of the total spin which 
are in principle observable in AC susceptibility or ESR 
measurements. Thermodynamic properties of AF sys- 
tems with uncompensated sublattice spins by now have 
been studied in great detail for vaijious anisotropy poten- 
tials and field configurations.ESo One main advantage 
of the small, high-symmetry modified FW studied in the 
present work is that the dependence of various thermo- 
dynamic quantities and spin correlation functions on the 
small number of microscopic parameters entering Eq. (|l|) 
can be evaluated analytically. 

The outline of this paper is as follows. In Sec. || we 
review the theory of spin tunneling in FWs and calcu- 
late the spin correlation functions Xaa- In Sec. |l| we 
discuss the modified FW within a classical vector model 
and show that, in contrast to the FW, tunneling of n now 



also results in oscillations of the total spin. In Sec. IV 



and Sec. |^, we develop a semiclassical theory of the mod- 
ified FW, thus generalizing earlier "warkjsn AF systems 
with uncompensated sublattice spinsEll Edto Bj, ^ 0. For 
I J'/ J— 1| ^ 1, we evaluate explicitly the tunnel splitting 
A and magnetization M^ of the modified FW, and show 
that J' can be determined by measuring M^- We calcu- 
late the Fourier transform of the real-time susceptibility 
x'zzi^ — ^/^) ^-nd prove that Neel vector tunnehng can 
be detected in AC susceptibility or ESR measurements. 
In Sec. [Vl[ we discuss the modified FW for the limiting 



present work for experiments is discussed in Sec. VII. 
Finally, we indicate that our results can be e asily g ener- 
alized to other systems such as ferritin (Sec. VIII ). We 



summarize our results in Sec. IX 



II. THERMODYNAMICS AND SPIN DYNAMICS 
OF THE FW 

A. Thermodynamics 

In this section, we give a brief review of previous work 
on Neel vector tunneling in FWs.Eil In particular, we 
point out that the notion of quantum tunneling applies 
only to Feio and Cs:Fe8 (and somewhat less to NarFcg) 
for large hx^ but not to Li: Fee. 

The minimal Hamiltonian of the FWs contains only 
two parameters, J and kz, 



Ho 



N 



N 



N 



'i+1 



+ hh-^s,-kzY^ s ■ 



(2) 



4=1 



4=1 



Here N — 10, 8, or 6 and s^v+i = Si, h = gfis^/K with 
B the external magnetic field, and g = 2 is the electron 
spin g-factor. Throughout this paper we restrict our at- 
tention to B = B^e^, i.e. magnetic fields applied in the 
ring plane. For kz = 0, the eigenstates of the total ^m 
S = J2i=i ^i are also energy eigenstates, with energyll3'0 
Es,s^ - {2J/N)S{S+i) + hhxSx. For systems with weak 
anisotropy, kz ^ 2J/{Ns)'^, the anisotropy can be taken 
into account in perturbation theory for a wide range of 
hx- However, the scenario changes for large anisotropy 
kz ^ 2 J/(iVs)^, where mixing of different spin multiplcts 
becomes appreciable, as is the case for Feio, CsrFcs, and 
NaiFee. 

Both the tunnel splitting A and the magnetization M^ 
of the FW can be obtained from the partition function 
Zo[hx] which we evaluate using spin path integrals. In- 
troducing sjjiij coherent states, we decompose the local 
spin fields s 



H 



(-ir 



1 



(3) 



into a Neel ordered field ±sn (n^ = 1) and fluc- 
tuations 1 _L n around it. For small systems con- 
taining 6, 8, or 10 spins, spatial fluctuations of the 
fields n and 1 are frozen out at low temperature T. 
Carrying out the Gaussian integral over 1, we obtain 

with a Euclidean La- 



Zo^JVn exp [- /„'^" dr Lo[n]/h 
grangean depending only on n. 



^o[n] 



~87 



[-(in X h - h)^ -I- (h • n) 



,2„2i 



(4) 



where ujq — sy/8Jkz/h. 



In contrast to the classical description of the states | t) 
and I I) used in Sec. [J, in a quantum mechanical treat- 
ment n always exhibits quantum fluctuations around its 
classical minima. The notion of quantum tunneling, how- 
ever, is only applicable if n is well localized in the states 
I t) and I i), i.e. if 1 - (t \nl\ ]) < 1. The quantum 
fluctuations of n can be estimated from Eq. (Q). For 
h = 0, Lq describes two independent harmonic oscilla- 
tors of frequency tJoj corresponding to fluctuations of n 
in the direction of Gx and e^. If the amplitude of the 
fluctuations is small, we can evaluate the mean devia- 
tion of the Neel vector from e^, 1 — (n^) ~ 2/{So/h), 
where Sq/H = Nsy^2kz/J is the classical tunnel action. 
Hence, n is well localized along e^ only if Sq/H ^ 2 or, 
equivalently, if the ground-state energy, 2 x wo/2, is small 
compared to the potential barrier NkzS^ . The scenario 
changes if a strong magnetic field h = /i^-e^, h^^ ojq, is 
applied in the ring plane. Then, the mode of n along Gx 
is frozen out, such that 1 — (n^) ~ l/{So/h), and the FW 
can exhibit quantum tunneling if So/h 3> 1. Note that 
for large tunnel action, So/fi ^ 10, the tunnel splitting 
becomes small, which would make the system under con- 
sideration a less favorable candidate for the observation 
of MQC. 

The FWs LiiFegp-NajEec—CsiFeg, and Feio have been 
well characterizedMlialla^El For Fcio, J = 15.56K and 
kz = 0|^88J. For CsiFcg, J = 22. 5K and k^ = 
0.019 1J.Ej For Fcg, J and fcz vary appreciably depending 
on the central alkali metal ion and ligands: for NaiFce, 
J = 32.77K and kz = O.OLS&i whereas for LiiFcg, 
J = 20.83K and k^ = 0.0053 JEHiJ In Table |, So/h and 
?iwo/2/is are given for Feio, CsiFcs, NaiFcg, and LiiFeg. 
As is obvious from these values, for none of the molecu- 
lar rings So/h is sufficiently large to assure that a tunnel 
scenario is rigorously apphcable if hx < wq. In Na:Fe6, 
even at large Bx ^ 20T, n is far less well localized along 
ez than in Feio, which hence remains the most favor- 
able candidate for the observation of quantum tunneling. 
Note however that even in Feio and CsrFeg, So/h is so 
small that corrections tOj-the instanton techniques used 
below may become large 
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So/h 


hLJo/2^lB 


Feio 


3.32 


7.68 T 


Cs:Fe8 


3.91 


16.37T 


NaiFee 


2.47 


20.11T 


Li:Fe6 


1.54 


7.98T 



TABLE I. So/h = Ns^/2k7JJ and hujo/2nB for Feio, 
CsiFeg, NaiFee, and LiiFee. 



For hx ^ wo: the magnetic field B strongly confines n 
to the (y, z)-plane and thus determines the tunnel path 
of the electron spins. This allows one to evaluate Zo- We 
parameterize 



(cos( 



, sm U cos ( 



(5) 



and expand Lq to second order in i!} = 6 — 7r/2, 
Lo[n] = ^-;-[-ihx ~ i<j>f " LolsiT? ( 



8J 



\dG-^W, 



(6) 

where G-i[0] = {Nh^ /^J){-d'^^ + {hx-i4>f ^ujlsix/^ (fy, 
4> = dr4>, and ujo = s\/8Jkz/h as defined above. The 
typical energy scales for the dynamics of cj) and -d are ujq 
and hx, respectively. Due to this separation of energy- 
scales, we can use an adiabatic approximation, in which. 
I? oscillates rapidly in a quasistatic harmonic potential 
{Nh'^/8J){h'l — 2ihx4>)- Integrating ouLj?, we obtain an 
expression for Lo depending only on 0.liJ For ksT /h ^ 
ujo ^hx, we find 



ioM 



Nh^ 
~8J" 



[^{hx^i^t^) — ^0 ^^'^ ' 
+O{iol/hx), 



i(j) 



(7) 



where the term h(hx — *</>)/2 arises from the § fluctu- 
ation determinant. The two saddle-points of Eq. (|^), 
(f) = 7r/2 and = 37r/2, correspond to the two classi- 
cal spin configurations sketched in Fig. ^. If tunneling 
is weak, So/h ^ ^nitije remaining path integral over (f) 
is straightforward.KtJ Summing all multi-instanton so- 
lutions, one finds 

'A^ft2 .- 



Zo = exp 



P 



-UJo 



8J 



cosh 



f P^jhx) 



with the tunnel splitting 
A{hx) = Ao 



sm 



Nh, 
n—h. 



(8) 



(9) 



where /3 — l/fc^T, Aq = 8hujo\/So/2Trh exp[— 5o/?i], and 
So/h — Nsy^2kz/J. In particular, A is periodic as a 
function of hx- Differentiating with respect to Bx, we 
obtain the magnetizationtj 



Mx = [giJLB) 



Nh 
4J 



1 aA 

2hdhx 



■ tanh 



/3A 



(10) 



Because dA/dhx is discontinuous at the zeroes of A, 
Mx exhibits steps at Bc^n — nAJ/Ngfis, where n = 
1,2,..., TVs. From 



Ma = (g^s)^^ / Vn [h^ - i(n x n),. 



n„h • nl 



xe 



-/.f drLo[n]/;i 



(11) 



it also follows that, for arbitrary hx. Ma = for ha — 0, 
a = y, z, which is a resrJ-t,of the invariance of Ho under 
rotation around B by ttEB 



B. Spin dynamics 

We show now that the tunnehng dynamics of n does 
not enter the susceptibihty of the FW, i.e. that the imag- 
inary part of the susceptibihty, Xaa{'^)i has no absorp- 
tion peak 6{uj — A//i). To prove this, we calculate the 

susceptibility of the total spin, S = X]i=i ^i' i'^ Matsub- 
ara representation, 

Xaa(*^„) = (5/ii3)' / dTe''^"^[(r,^„(r)^„(0)) 
Jo 

-{Sa.?l (12) 

where a — x,y, z. As we will show below, Xaa(j''^ri) con- 
tains no terms proportional to l/(iu;„ — A/fi). This im- 
plies that AC susceptibility or ESR measurements cannot 
be used to detect Neel vector tunneling in FWs described 
by ^0- 

Xxxii^n) can be calculated from the generating func- 
tional Z[6h.j;{T)], where Shx{T) is a small probing field 
added to the static field h^ ^ uiq. Because we are 
only interested in the low frequency response of the FW, 
LOn ^ A/h <^ ujQ, we may restrict our attention to a 
slowly varying field Shx whose Fourier components van- 
ish for ujn > A/h. The typical timescale for dynamics of 
(j), l/wo, is short compared to the timescale on which Shx 
varies. In particular, approximating Shxir) by a constant 
during instanton passage, we find 



-PE, 



^f3h 



['"" At 
Z[5hx{T)] ~ exp[y — 

hx + 6hx{T) +uj(i 



dr fNK- 



8J 



ihx+5hx{T)f (13) 



f^'^A^Aih^^^SKir)^ 



Differentiating twice and setting 5hx —^ 0, for tj„ < A/h, 
,[Nh 1 d^A 



Xxxiii^n) = (.gMs 

f3 fl dA 



AJ 2hdhl 

2 

-2 



tanh(/3A/2) 



ft V 2 dhx 



cosh-'(/3A/2)<5^„,o ■ (14) 



The transverse susceptibilities can be evaluated di- 
rectly from 5'^ Z / 5ha{T) 5ha{Q) , a — y,z. Using the 
parametrization in Eq. d^), i? — — 7r/2 can be inte- 
grated out in the path-integral expression for Xaa ■ After 
lengthy calculation (Appendix H), we obtain for aj„ <^ hx 



Xyy{i^n) 
Xzz{iuJn) 



t ^2^fi 

(ffA^s) ^> 



(15) 



From Eqs. (|l4| ) and (15) it is evident that none of the 
susceptibilities Xaaii^n) contains a term proportional to 
l/{iwn ± A/h). In the tunneling regime discussed here, 
\g) and |e) are energetically well separated from all other 
states. Then, Eq. (|l5|) and the spectral representation 



Xaa(iWn) ^^ ^^—\{i\Sa\i)\^ 



«J 



1 



1 



iuJn + {Ei-Ej)/h 



iujn - (El - Ej)/h 
-PhS^^ASc? (16) 



where \i) and \j) label energy eigenstates, imply that 
{e\Sa\g) = (a = x,y,z). Although for the parameters 
of Fcio, CsiFcg, and NaiFcg tunnel corrections to Xaa 
neglected in Eq. (09) can be significant, the main con- 
clusion of our calculation - that coherent tunneling of n 
does not enter the susceptibilities Xaa ~ remains valid. 

Indeed, {e\Sa\g) = is a direct consequence of the 
invariance of Hq as i ^ i + 1, i.e. the exchange of the two 
sublattices of the bipartite AF ring. In order to clarify 
this point, we introduce the sublatticc spin operators 



Sa = 2^ Si, Sb = 2^ 



(17) 



i— odd 



with S = S^ + Sb. In a semiclassical description of the 
FW, spins of one sublattice couple ierromagnetically to 
each other. The classical spin fieldsEJ s^ then obey s^ = 
Sa/{N/2) for odd i, and s^ = Sb/{N/2) for even i. 
amounts to treating S^ and S^ as single large spins 
with spin quantum number Ns/2, and Hq reduces to 



H{^> 



ih\ 



^±^ 

N^ 



Sb 
2k, 



f hh ■ {Sa 



Sb) 



^A,z 



^B,z) 



(18) 



Similar to the nonlinear sigma model (NLSM)-formalism 
used above, iJo,subi provides an appropriate description 
of the low-energy physics of the FW. 

In Fig.J|, the two classical ground-state spin configura- 
tions of i/o.subi are shown. A finite magnetic field B tilts 
the sublattice spins away from the direction of B such 
that S = Sa -I- Sb is parallel to B. During tunneling 
of n, the sublattice spins retain their position relative 
to each other and rotate jointly around e^,. The total 
spin vector S, however, remains invariant during tunnel- 
ing such that the real-time spin correlation function does 
not contain terms proportional to e /'' or, equivalently, 
(e|5„|5)=0. 





'B f A 

FIG. 2. Sublattice spins and total spin of the FW in mag 
netic fields. 



III. A PHENOMENOLOGICAL MODEL FOR 
THE MODIFIED FW 

The fact that the susceptibiUties Xaa of the FW show 
no signature of the tunnehng of n is a consequence of 
the symmetry of -ffo^suW under S^ ^^ Sb- In the modi- 
fied FWs [Eq. ™)] this symmetry is broken, which makes 
them much more suitable for the observation of tunnel- 
ing of n. In this section we discuss the saddle-point, 
i.e. classical, properties of the phenomenological sublat- 
tice pip^l for the modified FW. We generalize earlier 
workErEJ to finite magnetic fields and showj-that, in con- 
trast to systems with easy-plane anisotropy,E3 in molecu- 
lar rings with easy-axis anisotropy and finite excess spin, 
S oscillates as n tunnels. 

Following Sec. II B, we introduce S^ and Ss as defined 
in Eq. (0). The Hamihonian H of the modified FW 
can thus be mapped onto a simpler version in terms of 
sublattice spins, 

-Hsubl = 777 — , ^r N AT ^A • as + fih • (hA + SB) 



2kz 



{Ns + 2ds)Ns 
l + 2{kis''-k,s')/iNk,s^)^, 



(1 + 2Ss/Ns)^ 



S 



A.z 



s 



B.z 



(19) 



where Ss = s' - s, Sa = {N/2 - l)s + s' = Ns/2 + 5s, 
and Sb — Ns/2. Throughout this paper we assume the 
following inequalities: 



\5s\ <$:Ns, 

\k'^s'^ - k.s'^l <^Nk^s^/2,Js, 

2J's' < NJs, 



(20) 
(21) 
(22) 



where J, J' > 0. Eq. (GO) guarantees that the modified 
FW is an AF system with small excess spin, to which 
the theory of Neel vector tunneling applies. Both for 
Cr and Ga dopant ions {\ds\ — 1 and 5/2, respectively), 
Eq. (pfl) is well satisfied. Eq. ( pl| ) will allow us to treat 
the difference in sublattice anisotropics, [k'^s'"^ — kzs'^), 
in perturbation theory. Typical values of k'^ and k^ are 
on the order of only O.OIJ, such that this condition holds 
for most systems of interest. Finally, Eq. (|23) together 
with Eq. ( po| ) assures that the 'bulk' parameters of a FW 
are only slightly altered by exchanging one single spin, 
such that the parameters of the undoped FW [Table |] 
still determine whether the modified FW is in a quantum 
tunneling regime. However, as will be shown below, the 
excess spin Ss ^ leads to qualitative changes in both 
thermodynamic and dynamic quantities. We further as- 
sume 



fcsT ^ hujo, 



(23) 



which allows us to restrict our attention to the low-energy 
sector of the FW, which consists of two tunnel-split states 
only. 

We first discuss the classical vector model of Eq. ( p^ ) 
for kz — k'^ — 0, but finite h^- For an AF system with 



equal sublattice spins, the spins would lie close to the 
plane perpendicular to the field h^. As sketched in Fig. ||, 
tilting of the spins leads to a gain in energy Nh^h-^/8J. 
However, for uncompensated sublattices, the configura- 
tion sketched in Fig. pK provides an energy gain \6s\hhx 
and hence is favorable for hh^ <C \5s\ 8J/N — hhc- Only 
for hx 3> /ic, the classical ground-state spin configura- 
tion is as sketched in Fig. y|b. The energy is minimized 
if the projections of Sa and Sb onto the (y, 2)-plane are 
antiparallel, such that S = Sa + Sb is parallel to B. 
This picture remains valid for a system witi3-|easy-p/ane 
anisotropy and B applied in the easy plancEd 

The scenario changes for easy- axis (e^) anisotropy 
and a magnetic field B perpendicular to e^. First, 
the anisotropy favors the spin configuration sketched in 
Fig. Sb over that in Fig. pk. The true classical ground- 
state spin configuration depends on the ratio kz/hh^- 
More important, even for hx ^ he, S now has a com- 
ponent pcrpcndicuar to B. The reason for this is that 
for ^s 7^ or fc^s'^ — k^s"^ ^ 0, Eq. ( p^ is no longer 
invariant under exchange of S^ and Sb if kz 7^ 0. Due 
to Eq. (|2^) and (p2|) , the components Sy and Sz of the 
total spin can be evaluated perturbatively. With the po- 
lar angle parameterizing the projection of S^ onto the 



(?/, z)-plane, to leading order in (5s and (fc^; 



^)/2Js 



we obtain 



Sz 



''% 



kzs' 



2Js 



sm 



COS( 

sinf 



(24) 



As is evident from Eq. (M), finite Sy or Sz can be due 
to Ss ^ or k'^s'^ — kzs'^^0, i.e. excess spin or unequal 
effective anisotropics for S^ and S^. 

According to Eq. (pO|), the modified FW is an AF sys- 
tem with small excess spin which is expected to exhibit 
spin tunnel dynamics qualitatively similar to the FW, as 
indicated by the close formal analogy between Eqs. ( |lq ) 
and (pj|). For magnetic fields 



max[fiwo, \Ss\8J/N] < hhx < 4Js, 



(25) 



the sublattice spin vectors S^ and S^ lie close to the 
(y, z)-plane. Due to the easy-axis anisotropy, configu- 
rations with Sa and Sb close to the ^-axis are ener- 
getically favorable. It is noteworthy that the condition 
hx 3> Wo is not indispensible for quantum tunneling of 
n, but only assures that a tunnel sc enari o remains ap- 
plicable for a wider range of kz [Sec. II A |. In contrast, 
hhx ^ \Ss\8J/N guarantees that there are two energet- 
ically degenerate, macroscopically distinct spin configtu 
rations between which spin tunneling may take placecJ 
and hence, in general, will shift the range of magnetic 
fields in which a tunneling scenario as discussed in the 
present context [Fig. |l] is valid. Henceforth, we will al- 
ways assume that Bx is large enough to satisfy Eq. ([2^). 
For sufficiently large kz, the two-state model for the low 
energy sector of the system outlined in Sec. | then still 



applies. As Eq. ( J24| ) indicates, the modified FW exhibits 
one important novel feature: as n tunnels between e^ 
and — e^ (0 = 7r/2 and (f) — 37r/2, respectively), the 
z-componcnt of the total spin, Sz, oscillates between 
^0 = (Ssivi/hl + {k'zs'^ - kzS^)/2Js) and -So- 
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,-y 



■-y 



Sa S S^ 



(a) 



B 



B 



^, 



'-y 



(b) 



(c) 



FIG. 3. Classical ground-state spin configurations (sche- 
matically) of an AF system with excess spin in a magnetic 
field, (a) k^ = k'^ = and h^ < he ^ \5s\8J/Nh. (b) h^ > he 
and kz = k'^ = 0. (c) hx > he- Now, the finite anisotropics 
kz,k'^ 7^ lead to a finite component Sz perpendicular to the 
magnetic field B. 



IV. THERMODYNAMICS OF THE MODIFIED 
FW 

In order to quantify the statements on Neel vector tun- 
neling in modified FWs [Sec. [II |, we now develop a semi- 
classical theory of the modified FW. In this section, we 
discuss thermodynamic quantities such as the magnetiza- 
tion Mx and specific heat Cy. While our theory treats the 
spins semiclassically, we give the explicit dependence on 
the microscopic parameters of H [Eq. (||)] . To this end we 
evaluate the partition function Z , thus generalizing the 
procedure reviewed in Sec. [I A to systems with k'^ "^k^^ 
J' 7^ J, and (5s ^ 0. The most significant change is that, 
for bs 7^ 0, the staggering [Eq. (|3|)] must be modified in 
order to account for s^ = s'^ . The ansatz 



Sl 



s n 



Sj = (-l)*+isn + l V^T^l 



(26) 



is equivalent to the assumption that spins witiki. sjutlat- 
tices A and B are ferromagnetically coupled.l3'aEiro The 
results for the magnetization M,j. and susceptibilities ob- 
tained from this ansatz turn out to be in good agreement 
with those obtained from numerical exact diagonalization 
(ED) (see below and Sec. 0) as long as \J' j J - 1| < 1. 
We restrict ourselves to this case first and discuss further 
limiting cases J' f^ J and J' :^ J in Sec. VI. 

As for the FWjlJ at low temperatures spatial variations 
of the Neel field n, and the fluctuations 1 around it, are 
suppressed in small ring system. With the coherent spin 
states defined in Eq. (E§),L3 



Z^ fvnV\S{nl)e~ 



S[n,\]/h 



(27) 



where 



S[n,l]^-ih[s'u;[n- 



ill 

s 



N 



^E^[(-i 




l3h 



drHlnA] 



(28) 



Here, 



H[n,l] = 2[{N-2)J 
-[{N-I)kzs^ 



2-J']l2 
s 



h{N 



k'^s'^]72l + -{kzs^ - k'^s'^)nJz 



-Ss hh ■ n 



(29) 



is the classical energy of a given spin configuration. The 
small term [{N — l)kz + k'z]l1 has already been neglected. 
The first term in Eq. (g8|) is the Berry-phase term, where 
w[ri] denotes the area traced out by some vector ^{t) on 
the unit sphere. Here, we use the south pole gauge (i.e. 
uj[^{t) = —Gx] — 0). Expanding the Berry-phase term 
to leading order in 1 with the parametrization Eq. (0) 
yields 



I3h 



dr 



i5s(j){l -t-cos6') +i{N A )(n x n) • 1 



(30) 

The first term is due to the fact that, for 5s ^ Q, 
the Berry-phase terms of the AF ordered components 
s'w[n] + sS,^2'^[(-l)*+in] = /dr(7Vs(i4^0(l-Hcos0)) 
do not add to an integer multiple of 27r.EZlH 

Carrying out the Gaussian integral over 1 we obtain 
Z — Jl?nexp[— J„ dTL[n]/fi,], with a Euclidean La- 
grangean 

rr 1 N^^ , n , • •n2 

L[n\ = — ^[— (h + hATizez — in x n) 



8J 



where 



+ ((h + hArizez) ■ nf - ui^nl 
+Ssh[h.- n-i(j){l + cos 61)], 



J__ {N-2)J + 2s'J'/s 

N ^ 



(31) 
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OJ, 



(N + 5s/s)^ 
2{kzs^ - k'^s'^) 
{N + Ss/s)sh ' 
8J, 

'n 



(32) 



[{N ^ l)kzs' + k'. 



Eq. (|l|) is the analog of Eq. (|) for the FW. A finite ex- 
cess spin Ss of the modified FW leads to two significant 
changes in L[n]: first, the typical energy scales of the 
system are slightly renormalized, J ^ J and ljq ^ uj^^ 
even for J' = J and k'^ = kz- More importantly, L 
acquires an additional term due to the Zeeman energy 



and Berry phase of the uncompensated spin Ss n. For all 
cases of experimental interest, due to Eq. (|2l| ) we have 
hA <S^ oJo ^ hx- The ft^-dependent terms in Eq. 1^^) 
hence lead only to minor modifications of the thermody- 
namic properties of the modified FW compared to the 
FW, but feature in the dynamics. 

We will show now that, for h^ 7^ 0, n and hence the 
excess spin no longer trace a tunneling path in the (y, z)- 
plane. As can be seen from 



u 



^^[^{hx ^ 'i'P) — ^0 ^i^^ (l)]—ihSsi 



8J 



8J 

52 I „2 



+ cos 9 [{hx - i(p) + uJq sin"^ 



4J 



+2 cos 6 (hx — i4>) I 6s — ^ + Ha sin^ 9 sit? i 
\ Nn 



(33) 



the timescales for the dynamics of cj) and 9 separate due 
to Eq. (p5|), and we can again invoke the adiabatic ap- 

9 oscillates in a slowly 



[IB 



proximation used in Sec. 

varying harmonic potential with the potential minimum 

at 9o, where 



cos 00 — — 



hx - iij) 



{hx - i4>Y + 



a>n sin 



4^ . ^ . 2 
— ds + Ha sm i 



(34) 



Corrections to the adiabatic approximation are beyond 
the order of the present calculation. Eq. (34) shows that 
finite 5s ^ Q or Ha ^ Q leads to a shift in the saddle- 
points of the Lagrangean L[(l),6] away from 9q = 7r/2, 
which is due to the fact that then Sa — Sb no longer lies 
in the {y, z)-plane [Fig. He]. Expanding Eq. ( p3| ) to second 
order in i? = 9—9q and carrying out the Gaussian integral 
over 1!}, we obtain a 0-dependent effective Lagrangean 

L[0] = — —[—{hx — iij))'^ — u}q sin^ (j)] — ihlSs + - 
SJ V ^ 



Using the same techniques as for the FW, we find 



Z = exp 



/3('^/.^-;.^ + '^" 



V 8J " ■^ 2 
with the tunnel splitting 



coshf^J, (36) 



A{hx) = Ao 



i — ;,„ 4- Ss 



smTT I — -hx 
V4J 



(37) 



where S/h = NLbo/2J, Aq = Shujod S /2'Khexp[—S /h]. 
From Eq. (^^ it is also straightforward to derive all ther- 
modynamic quantities of interest. In particular, for the 
free energy F, the magnetization Mx, and the specific 
heat Cy we obtain 



8J 



1 



F = h " : •"" - ^-^^hl ~ ^ In cosh 1^1, (38) 



Mx = (.gMs) 



Nhh 



1 1 9A , I (3 A 

— = ! — — tanh 

4 J 2 2hdhx \ 2 



cv = ks i -^ cosh ^ — 



(39) 
(40) 



The most significant change in the thermodynamic prop- 
erties of the modified FW is that, for half-integer Ss, 
the zeros of A and hence the magnetization steps are 
shifted by a magnetic field 2J/NgfiB, i-e. half of a mag- 
netization plateau, compared to the unmodified wheel 
[Fig. Q. The magnetization plateaus then lie at half- 
integer spin values. At low T, the specific heat Cy should 
exhibit a characteristic Schottky anomaly, with a peak 
at Tq ~ 0.4A/A:b- So far, c„ has been measured only for 
Bx — in Fcio and Fee samples.l23 Measurements of Tq 
in Cy (T) for various hx ^ ojq would in principle allow one 
to observe the characteristic sinusoidal variation of A as 
function of hx- 



^h^+0{Cjl/hx,NhACja/8J). 



(35) 



Comparison with the corresponding Lagrangean for Ss — 
[Eq. (1^)] shows that, to leading order in the excess 
spin Ss and anisotropy field hA, the only effect of an ex- 
cess spin is to introduce an additional topological phase 
—iSscj). In particular, in contrast to the case hx = 0-4^^ 
cussed in earlier work on tunneling in ferrimagnets,LlrEj 
the potential barrier and hence the real part of the tun- 
nel action is only slightly altered by the excess spin. This 
is due to the factEa that, for hhx > |^s|8J/iV [Eq. (H)], 
the system is in the AF regime in which the tunnel split- 
ting is only slightly modified by the excess spin. Note 
that Eq. (35) is formally identical to Eq. (0), which pro- 
vides a rigorous proof of the statement that the modified 
FW also may exhibit tunneling of n for sufficiently large 
anisotropy, as already claimed on basis of physical argu- 
ments at the end of Sec. Ill . 
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FIG. 4. Schematic plot of the ground-state magnetization 
of a modified FW with integer ( — ) and half-integer (- -) 5s, 
respectively. 

The width of the magnetization plateaus, SBc,n — 
Bc,n+i — Bc^n = "iJ /NgiiB, whcre Bc^n is the field at 



which the magnetization M^ exhibits the n step, is 
a quantity which is accessible in experiments and from 
which J' can be inferred. Our theory predicts 



^Bc,„(J')-'5S,,„(J'- J) = 8 



s'ls {J' - J) 



[N + 6s/ s)^ gfiB 



(41) 
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In Fig. Bwe compare the functional dependence predicted 
by Eq. pl|) with the results of exact diagonalization (ED) 
on small rings, N = 4, s = 2, s' = 3/2 (upper panel) aij 
s' = 1 (lower panel), respectively, for \J'/J — 1| < 0. 
The ED results are in good agreement with the analytical 
result. The deviations for | J'/J — 1| > 0.1 signal the 
breakdown of our ansatz in Eq. (p6|) for J' significantly 
different from J (see Sec. VI below). 
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FIG. 5. Comparison of ED (symbols) and analyti- 
cal ( — ) results for the difference in plateau width, 
5Bc,n{J') — 5Bc,n(J' ~ J) as function of J' — J for Ss = —1/2 
(upper panel) and 5s = —1 (lower panel). A'^ = 4, s = 2, 
kz = k'z = O.IJ. The numerical error of the data points 
±0.002 X AJ/NgnB is smaller than the symbol size. 

We show now that although, for a given direction 
of n, the total spin S acquires a component perpen- 
dicular to the field B, the magnetization Ma = 
still vanishes for a = y, z. We define the fields 
inair) = S J dTL[n]/5{hha{T))\h^=o such that Ma = 
-{gfiB) JVnma{T)e~ f '^''^/'''/Z. Using Eq. (|3|), we ob- 
tain 

7712 = —j[-hAnz{l - nj 

+hxnxnz + i{n x n)z] + 5s Uz 



-riASiTKp + ds 

4J 



{h^ - i(t>y 



Nh ■ 
4J 



(42) 



with do defined in Eq. ( |34| ) .EJ As follows from the invari- 
ance of H under rotation around B by tt. My — M^ = 



for arbitrary h^- In particular, at T = 0, M^ — in- 
dicates that the ground-state is not a state with definite 
direction of n, but rather a coherent superposition of such 
states, (I T) + I i))/V2, as expected for a system which 
shows coherent Neel vector tunneling. 



V. DYNAMICS OF THE MODIFIED FW 

As we have shown above, the effective action L[ip] of 
the modified FW [Eq. (|35|)] is formally identical to that of 
the FW [Eq. (R)]. In particular, for large anisotropy k^, 
such that S/h ^ 1, th e mo dified FW is in the quantum 
tunneling regime [Sec. II A]. In this section, we evaluate 
explicitly the spin susceptibility X2z(t) for the modified 
FW. 

In order to motivate this, we first calculate Xzz using 
the results of the classical vector model [Sec. III]. For n 
along Bz or — e^, the z-component of the total spin vec- 



tor is finite, Sz — ±5'o, where S'o = Ssujo/hl + {k'^s'^ 
kzs'^)/2Js. For n(t = 0) = e^, the coherent tunneling 
of n then results in an oscillating Sz{t) — S'o cos ( At/ /i), 
such that the Fourier transform of the (real-time) sus- 
ceptibility exhibits an absorption pole x'zzi^ — ^/^) = 
Tr\So\'^ tanh{pA/2)S{uj - A/h). 

Generalizing the procedure for the FW [Appendix H] 
to (Ss 7^ 0, |WE calculate the quantum corrections to this 
result fromE3 

Xzz{t) = (5^b)2^(1 - (sin2 0))5(T) 

with mzir) given in Eq. (0). As for the undoped FW, 
the correlations of the ^^-terms in mz^r) give rise to 
a stronglyp4)eaked term (g/iB)^iVM(T)(sin" 0)/4J. We 
hence fincO 



Xzz(T) = (g^J.B) TT'^M 



4J 



(44) 
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In stark contrast to the FW, the path integral in Eq. (|4J) 
gives rise to terms proportional to exp[±A|T|/fi] such 
that, upon Fourier transform, the susceptibility in Mat- 
subara representation contains terms l/(ztt'„±A/fi,). The 
path integral is most easily evaluated in a Hamiltonian 
description. We requantize the field (f) and use an effective 
two-state Hamiltonian to evaluate the matrix elements. 
Inserting the expression for Ha, we find [Appendix y[ 



(elSz 



N{k'^s'^ - kzs^) 
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2Js{N + ds/s) 
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-sSU 



J 



hi y "NhCjQ 
x'L{^'^K)^^{g^,Bf\{e\S^g)\^ 

'pa 



X tanh 



5{oj~A/h). 



(45) 



(46) 



Eqs. (^a) and (M) are the main results of this section. 
For 6s = —5/2 or —1 (for Ga and Cr dopants, respec- 
tively), KejS'zjg)! can be of order 0.1 even for hx > Swq- 
For ksT < A the susceptibility of the modified FW then 
exhibits a resonance at oj = ztA/h which is accessible 
in AC susceptibility or ESR measurements. The terms 
J/NhhjQ = l/{2S/h) in |(e|5z|5)| are quantum correc- 
tions to the classical result derived at the beginning of 
this section. 

So far we have ignored decoherence of the spin tunnel- 
ing, which is crucial for the notion of MQC. The condi- 
tion r ~ A/fi,, where T is the electron spin decoherence 
rate, marks the transition from coherent to incoherent 
tunneling dynamics. As is evident from the classical vec- 
tor model discussed in Sec. II], Sz follows the tunneling 



dynamics of n. In particular, for a single modified FW, 
the decay rate of \{nz{t)nz)\, F, is also the decay rate of 
\{Sz{t)Sz)\. For F ^ 0, the (5-peak in Eq. (H) is then 
broadened into a Lorentzian of width F. In experiments 
carried out on an ensemble of modified FWs, inhomoge- 
neous broadening (e.g. due to crystal defects or nuclear 
spins) adds to the width of the resonance peaks. The ex- 
perimentally determined linewidth of the absorption and 
emission peaks provides an upper limit for F. This should 
allow one to settle the experimentally unresolved prob- 
lem of whether true quantum coherence is established in 
ferric wheels. 

Finally we compare our result for the transition matrix 
element |(e| 5*2 1.9)1 entering Eq. (^6|) with results obtained 
from ED on small rings for a wide range of parameters 
[Figs. |g, 0, ||. For simplicity, we assume J' = J. In 
the range of validity of our theory [Eq. (p5|)], for 6s ^ 
[Figs. ^, |g], the agreement of ED (o) with analytic re- 
sults (-pi is both qualitatively and quantitatively con- 
vincing.Cj The small oscillating features seen in the exact 
results are due to tunnehng corrections 0(exp\—S/h]) to 
Kel^zlg)!, which were neglected in Eq. (|45|). For Fig. 0, 
where 6s = 0, our theory makes the correct qualitative 
prediction that KelSzl^)] depends only weakly on h^, but 
overestimates the matrix element. However, due to the 
smallness of the matrix element for 6s — 0, the discrep- 
ancy can be due to terms neglected in the derivation of 
Eq. (pq). The significantly different qualitative features 
of FigT^ compared to Figs, g and ^ arise from the fact 
that ^s = in Fig. ^ The difi'erent functional depen- 
dence of KejS'zIg)! on hx for the two cases 6s ^ and 
(5s = is well understood within the theoretical frame- 
work presented here [Eq. (|45|)]. The very large difference 
in matrix element magnitude illustrates the importance 
of looking at doped rings in experiment. 



We conclude this section by remarking that, for fi- 
nite excess spin 6s, the second transverse susceptibility 
Xyy{^) also has an absorption pole a,t uj — A/h. How- 
ever, since Gj, is a hard axis, the spectral weight of this 
pole is significantly smaller than that ofx'izi^)^ such that 
Neel vector tunneling in the modified FWs can be more 
easily detected by probing the latter quantity. 
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FIG. 6. Transition matrix element | {e\Sz\g) \ for small rings, 
TV = 4 with J' = J, s = 5/2, s' — 2. In the upper panel, 
kz — k'z = 0.1. In the lower panel, kz — k'^ = 0.055J is 
chosen such that S/h ~ 3.3 as for Feio. The analytical result 
( — ) is shown for max[ficio, |(5s|8J/A'"] < hh^ < 4Js. Due to 
Eq. (Eq), our theory is rigorously valid only for fields much 
larger than max[hijJo, \Ss\SJ/N] and much smaller than 4Js. 
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FIG. 7. Transition matrix element KelSzIg)! for a small 
ring iV = 4 with J' = J, s = s' = 2, i.e. 5s = 0, but 
k'z — 1.5kz and hence Ha j^ 0. Again, kz — 0.0655J is chosen 
such that S/h ~ 3.3, as for Feio. The analytical result ( — ) is 
shown for inax.[hijJo, \Ss\8J/N] < hhx < 4Js. See also caption 
of Fig. I 
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FIG. 8. Transition matrix element |(e|52|g)| obtained with 
the phenomenological sublattice Hamiltonian [Eq. (119)] for 
Feio with one s — 5/2 substituted by (a) a dopant with 
s' = 3/2 (e.g. Cr) (numerical data o, analytical prediction 
— ), and (b) a dopant with s' — (e.g. Ga) (numerical 
data o, analytical prediction - -). Note that, in this case, 
the numerical data is not obtained from ED of Eq. (hi), but 
rather of Eq. (|l9|). For simplicity, we assumed J' = J and 
k'^ = kz — 0.0088J. The analytical results ( — and - -) are 
shown for max [too, |5s|8J/A''] < hhx < 4Js. See also caption 
of Fig. |. 



VI. THERMODYNAMICS AND SPIN 
DYNAMICS FOR J' /J > 1 AND J' /J < 1 

The deviations of the ED results from our theoretical 
predictions shown in Fig. ^ indicate that^for J' /J ^ 1 
or J' /J ^ 1, the theory developed in Sec. |^ is no longer 
immediately applicable. Indeed, results obtained by ED 
for the ground-state magnetization M^ in small rings 
(iV ^ 4, s = 2, and s' = 3/2) [Fig. || indicate that 
one of the main results of Sec. |^, that M^ exhibits a se- 
ries of equally spaced magnetization steps, does not hold 
any more. As we will show below, this is due to the fact 
that our ansatz Eq. ( |26| ) needs to be modified for J' sig- 
nificantly different from J. In this section we show that, 
for the limiting cases of J' 3> J or J' ^ J, the modified 
FW can be mapped onto the problem discussed in the 
preceding sections. We discuss the qualitative features 
of Mj, for these systems and show that coherent tunnel- 
ing of n also results in coherent oscillations of the total 
spin. 

J' ^ J : In this limit, Si decouples from all other spins 
and aligns antiparallel to B for hhx ^ J's. The remaining 
spins S2, S3, . . .sn form an open spin chain, as sketched 
in Fig. Il^a. As shown in Appendix O, the Lagrangean 
of an open spin chain with an odd number of spins can 
also be mapped onto Eq. (Esf), with Ss = s and slightly 
renormalized J = JN{N - 2)/{N - l)^. We predict that 
Mx has the following features: 

• Mx > gfJ-sis' + s) for 2 J's < hh^. 

• For Taa,x[huia, s8J /N] <C hhx "^ 4Js, Mx exhibits 
a series of equally spaced magnetization steps with 
plateau width SBc^n — 'iJ/Ng^B- Depending on 



whether 6s is half-integer or integer, the plateaus 
correspond to states with half-integer or integer to- 
tal spin, respectively. 

J' ^ J : In this limit, the spins s^v, si and S2 are 
strongly coupled. In a semiclassical picture, Si aligns 
antiparallel to sjv and S2 and the three spins act as one 
single spin |2s— sjcoupled to S3 and Sjv-2 with exchange 



constant J [Fig. |10|b]. For simplicity we assume s' < 2s. 
Then, for hhx ^ J'{s — Ss + 1), H can be mapped onto a 
Hamiltonian of the form Eq. (|I|) with N -^ N-2, J' ^ J, 
k'^ -^ (fc^s'2-h2fc^s2)/(2s-s') , and s' -^ 2s-s' = s-6s. 
Because all iV — 2 exchange couplings in the new Hamil- 
tonian are identical, the theory developed in Sec. IV and 



|y| remains applicable. In particular, for the ground-state 
magnetization Mx we make the following predictions: 

• For ma,x[hCjo, \6s\8J/{N - 2)] < hhx < 4Js, Mx 
exhibits a series of equally spaced magnetization 
steps with SBc^n — 4J/(iV — 2)g^B- Depending on 
whether 5s is half-integer or integer, the plateaus 
correspond to states with half-integer or integer to- 
tal spin, respectively. 

• For hhx ^ J'{s — 5s + 1), the Zeeman energy is 
sufficiently large to destroy the AF configuration of 
sjv, Si, and S2. This results in a series of additional 
magnetization steps with spacing J' . 

Note that a similar argument also applies if J' < 0. In 
this case, the three spins sjv, Si, and S2 are ferromagnet- 
ically coupled and align parallel. Again, the system can 
be mapped onto a smaller ring (as in Fig. |lO|b), where 
now s' — > 2s + s'. The magnetization curve resembles 
the one shown in the upper panel of Fig. ||. In Fig. ||, 
ED results for small rings with N — A, s — 2, s' — 3/2 
are displayed. The qualitative features agree with all the 
above predictions. 
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FIG. 9. M^{B^) at T = for a small system with J' <^ J 
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FIG. 10. (a) For J' <^ J, si decouples from all other spins 
and the modified FW can be mapped onto an open spin chain 
excluding si. (b) For J' ^ J, sjv, si, and S2 are strongly 
coupled such that they can be described as one single large 
spin. The Hamiltonian of the ring then maps onto that of a 
modified FW with J' = J. 

We conclude that even the qualitative features of M^ 
allow one to estimate the parameter J' of a modified FW. 
Even more important, as we have shown, also for J' ^ J 
and J' ^ J the modified FW can be mapped onto the La- 
grangean [Eq. (31)] of a system which exhibits quantum 
tunneling of n. In all cases discussed above, the quantum 
tunneling of n leads to coherent oscillations of the total 
spin S, and thus can be observed in AC susceptibility or 
ESR measurements. 



VII. DISCUSSION 



The theory described in Sees. |l^, |^, and VI allowed 
us to derive explicit expressions for both thermodynamic 
quantities [Eqs. (|39| ) and (^] and the susceptibilitiy x'L 
[Eqs. (Esl) and (|46|)] of modified FWs. In order to estab- 
lish a connection with experimental issues, we now out- 
line the steps necessary to detect coherent Neel vector 
tunneling. For simplicity we restrict our considerations 
to Feio with Ga {Ss = —5/2) or Cr {5s = — 1) impurity 
ions, and assume J' ~ J. 

For finite excess spin (5s, the two energetically degen- 
erate spin configurations [Fig. |l| required for coherent 
spin tunneling as discussed in the present work certainly 
exist if hh^ > \Ss\8J/N [Eq. (|25|)]. This tunneling 
regime is well within experimental reach for Cr dogants 
{B^ > 9T), but not for Ga dopants {B^ > 23T).H For 
Cr dopants (s' ^ 0), however, the two new parameters J' 
and k'^ introduced in H [Eq. ffl)] must first be determined 
in order to characterize the ring system. 

Both J' and k'^ can be obtained from the measure- 
ment of two independent thermodynamic quantities, such 
as the ground-state magnetization and tunnel splitting. 
A schematic plot of the ground-state magnetization for 
integer 5s (o) is shown in Fig. 0. Although the mag- 



netization steps are smeared out at finite temperature, 
for T < IK, the magnetization step spacing 5Bc^n still 
can be obtained with high accuracy.ll3 With (5i?c,n = 
4:J /NguB and Eq. (^), this allows one to determine 
J'. The on-site anisotropy A:^ can be obtained from 
A, which depends sensitively on the tunnel action S oc 
\/{N — \)kzs'^ -t- fc^s'2 and hence on k'^ . The tunnel split- 
ting A (and hence k'^) is accessible either in AC suscep- 
tibility or ESR measurements [Eq. ([46|)], or in measure- 
ments of thermodynamic quantities, such as c^. Torque 
magnetometry is another experimental technique which 
has been used to determine tbe, anisotropy constant with 
quasi-spectroscopic accuracy.Ell 

Once J' and /c^ are known, Eq. (Hfl) determines both 
the position and the spectral weight of the resonance 
in x'zz{^) which arises from coherent quantum tunnel- 
ing of n. The characteristic functional dependence of A 
[Eq. d?!)] and \{e\Sz\g)\'^ [Eq. (H)] on K predicted by 
our theory can be checked experimentally. Finally, it is 
noteworthy that, although A can be determined from 
thermodynamic quantities, the key problem of MQC is 
the measurement of the decoherence rate F which is ac- 
cessible only in dynamic quantities, such as AC suscep- 
tibilities. 

Throughout the current work, we have considered FWs 
with only one dopant ion. As we have shown in the pre- 
ceding sections, thermodynamic and dynamic quantities 
of doped FWs may differ significantly from those of un- 
doped FWs. In the large samples investigated experi- 
mentally, doping will lead to a random distribution both 
of the number of dopant ions and of their position rel- 
ative to the direction of the magnetic field. We defer a 
detailed analysis of these issues to a future publication. 
Here we note only that the random distribution of impu- 
rities does not invalidate the above considerations, and 
stress the qualitative features which ensure this. The 
choice of a low impurity concentration results in a large 
majority of the FWs containing no dopants or having 
only one dopant ion, which allows one to extract the sys- 
tem parameters of the singly doped FWs. When intraring 
dipolar interactions make a significant contribution to the 
effective uniaxial anisotropy k^ , doping with only one ion 
changes the effective anisotropy from uniaxial to biaxial. 
The theoretical framework presented in this paper can 
be readily extended to account for biaxial anisotropies. 
Because the original uniaxial anisotropy dominates the 
biaxial correction, the altered tunnel splittings A^ in a 
singly doped FW have a magnitude similar to A for the 
undoped FW, and a separation which is small by compar- 
ison. Thus AC susceptibility or ESR measurements can 
be expected to observe signals corresponding to reversals 
of the total spin accompanying Neel vector tunneling, 
and governed by the frequencies A^. 
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VIII. SPIN QUANTUM TUNNELING IN 
FERRITIN 

The theoretical framework developed in this work is 
quite general and applies to other systems besides AF 
ring systems. In particular^ the results of the classi- 



cal sublattice model [Sec. [11 1 can be easily extended 
to different systems. In order to illustrate this point, 
we now discuss natural horse-spleen ferritin and artifi- 
cial magnetoferritin, in which spin quaBtum tunneling 
has alre£u±y_been studied experimentallytS c2l and theo- 
retically.c3c3c3 The experiments were carried out in the 
presence of small static magnetic fields (B^ < 10~^T). 
In this regime, S^, ^ leads to an energy bias between 
the states | t) and | J,), and tmmcling is suppressed for 
increasing B^. r— ■ 

Recently it was showntJ that, for sufficiently large 
field, i.e. hh^ ~ \Ss\ 8J/N, there are again two energeti- 
cally degenerate spin configurations between which spin 
quantum tunneling may take place. In natural horse- 
spleen ferritin^ J - 200K and Ss/N ~ 0.05. For 
a system with uniaxial hard axis anisotropy, Han.z = 
kz X]i sf z7 the tunnel barrier and hence the tunnel action 
can be effectively controlled over a wide range of parame- 
ters by-j/arying the magnetic field B^ applied in the easy 
plane." Tunneling in the plane perpendicular to BLgi: 
rise to a topological phase acquired by the spins .eZIc2 
A drawback of the setup considered in Ref. ^ is that, 
for uniaxial hard axis anisotropy, spin tunneling leaves 
the total spin S invariant if B is applied perpendicular 
to the hard axis. As for the FW with equal sublattice 
spins, spin tunneling cannot be observed in AC suscepti- 
bility measurements. 

However, experiments indicate that, in addition to 
the strong hard-axis anisotropy, ferritin also exhibits a 
second weak hard-axis anisotropy Han.y — ky^^sf , 
wherec3 ky/kz ~ 10^'^. In self-sustaining films of natural 
horse-spleen ferritin, the hard axis e^ is perpendicular to 
the film.E3 In the simplest experimental setup, interfer- 
ence of different spin tunnel paths then could be explored 
if a field B = B^e^ is applied along the hard axis. As 
long as hhz ^ NkzS^ /5s (~ lOK for horse-spleen fer- 
ritin), due to the large anisotropy energy the spins re- 
main confined to the film plane such that there are again 
two energetically degenerate spin configurations similar 
to Fig. n^. Tunneling takes place in the plane perpendic- 
ular to B, with a tunnel splitting 



A = Ao 



COSTT S\ 



Nh 

u 



(47) 



where Aq — 5 x 10~^K and the total staggered spin is 
S'tot — 2.5 X 4500 for natural horse-spleen ferritin. A 
is periodic as a function of Bz, with a period SBz = 
AJ/NgfXB (— 0.13T for natural horse-spleen ferritin). 
The advantage of this tunnel scenario is that quantum 
tunneling of n also results in a tunneling of the excess 



spin (5s, and hence leads to a large resonance peak in the 
susceptibility per ferritin molecule, i.e. 

X'L{^ c, A/h) c, ^n{Ssr tanh (^^) 5{u; - A/h). 

(48) 

The factor 1/2 takes into account the random distribu- 
tion of easy axes in the film plane. Due to the spread in 
particle number, the total staggered spin of the system 
can be either integer or half-integer. Hence, one will ob- 
serve two different tunnel splittings varying with hz as 
A = Ao|cos(7rA^;i/i2/4J)| and A = Ao| sin(7rA^;i/i2/4J)|. 
An experimental confirmation of this behavior would pro- 
vide further strong evidence that the resonaiiee observed 
in AC susceptibility measurements in ferritirJlSc^ is due 
to macroscopic spin tunneling. In particular, the period 
of the oscillations of A as function of the applied field Bz 
would allow one to estimate the total number of tunnel- 
ing spins. 



IX. CONCLUSION 

The AF ring systems discussed here, modified ferric 
wheels which are already available to experimentalists, 
combine the advantages of AF and ferromagnetic molec- 
ular magnets. The tunnel splitting A is sufficiently large 
that quantum coherence between macroscopically differ- 
ent states is established. Tunneling of the Neel vector n 
also leads to a tunneling of the total spin S, thus making 
the spin dynamics in modified FWs accessible to exper- 
iment. We have considered the simplest realistic model 
Hamiltonian H [Eq. (|])] for a system in which, for ex- 
ample, one of the Fe ions is exchanged by Cr or Ga. We 
showed that the additional parameters entering H can 
be inferred from equilibrium quantities such as the mag- 
netization. Moreover, for a wide range of parameters, 
the system still exhibits MQC in the form of coherent 
tunneling of n. Finally, we calculated spin correlation 
functions of the modified FW and showed that tunnel- 
ing of n can indeed be observed in AC susceptibility or 
ESR experiments, which allow one to measure both the 
tunnel splitting A and an upper bound for the spin de- 
coherence rate F. Hence they should be appropriate to 
verify experimentally that Neel vector tunneling in FWs 
is coherent. 

Throughout this work we have used spin coherent- 
state path integrals leading to a NLSM-description for 
the modified FW. The main advantage of this technique 
over ED is that thermodynamic and dynamic quantities 
can be evaluated for a realistic system size. In addition, 
an intuitive physical understanding of the spin dynamics 
(quantum tunneling of n) can be obtained. A drawback 
of the analytical approach chosen in the present work is 
that it naturally requires approximations. Corrections 
to our results, in particular 1/s-corrections, may become 
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appreciable for the parameters of the FWs. However, our 
analytical results for Kel^zl^)] agree well with ED results 
obtained for small systems. For the parameter range ex- 
plored in ED, deviations from our theoretical predictions 
become significant mainly if s' is small (s' — 1/2 or 1), 
where the ansatz Eq. (|2^) fails due to large quantum 
fluctuations of §1. Although numerical work on rings 
with iV = 6, 8, and 10 is challenging, some results have 
been obtained for the ground-state magnetization and 
torque.Ej A detailed numerical study also of spin corre- 
lation functions would provide clearer evidence for the 
range of validity of the present approach and could ex- 
plore its limitations. For example, as is well known, the 
sublattice Hamiltonian Eq. (jlj) is-epcact for A^ = 4. For 
increasing A^, however, ED resultscj show that there are 
small deviations from -ffo,subi- 

Note that our work also has importantJmplications for 
undoped ferric wheels. |_B,ecent torque^3 c^,E2I and pro- 
ton l/Ti-measurementstiJon single crystals of various Feg 
compounds indicate that these FWs could exhibit physics 
beyond the Hamiltonian Eq. (0). One important future 
step in explaining the new experimental data will be to 
clarify to which extent the observed phenomena can be 
attributed to inhomogeneous level broadening. The the- 
oretical framework presented here allows one to calculate 
analytically the inhomogeneous level broadening result- 
ing from a random distribution of single exchange cou- 
plings J' and on-site anisotropics k'^ which could be a 
consequence of lattice defects in Fee crystals. Inde 
cent work on Mni2 suggests that lattice distortiong|£2 j 
a distribution of anisotropy energies and (7-factorst 
could account for the observed broad distribution of tun- 
neling rates in Mni2. 

Finally, we stress once more that the Hamiltonian H 
[Eq. (|l|)] discussed in this paper is a simple model Hamil- 
tonian, which still leads to fascinating novel features in 
the physical proper ties of the modified FWs. However, 
as discussed in Sec. VH, realistic systems might require 
modification of Eq. (y). Generalization of the present 
approach to more complicated anisotropics is, however, 
straightforward. 
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APPENDIX A: SUSCEPTIBILITIES 

In this appendix, we sketch how the transverse suscep- 
tibilities Xaai'ii^n), a = y, z, Can be evaluated for the EW 
with Hamihonian Hq [Eq. @]. With My = M,_ ^ 0,0 



-(5Mi 



I Vnl^in.n). +n.h.n.]. 



Nhy 1 

[i{n X h)„ +n„;^,n,]oe-^'"'^^^°["l/^ (Al) 



where the first (second) square bracket is evaluated at r 
(0) . Using the parameterization in Eq. (m and expanding 
to second order in 1!) = 6 — n/2, for a = y the square 
bracket reads 



i(n X h)y + UyhxTi^ 



-{hx — i4>)'d cos (j) — id sin 0. 

(A2) 



The corresponding expression for a — z can be obtained 
by setting — > — 7r/2. Integrating i9, we obtain 

Xyy{r) = (5Mb)'^(1 " {cos" 4>))6{t) 



V(j)e-^o^^-Lol4>]/n 



Nhy 1 
4jy z 

X [G^^{T){hx - i0(r)) cos ^(r) (/i^, - iq 
+iG^^{T){hx -iij){T)) cos (f){T) sin(^ 
+iG^^ (r) sin 0(t) {hx ~ icj)) cos </) 
+*^^^,j ('^) sin '/'(•r) sin 0] . 



I cost 



(A3) 



The Green functions are defined by G§^{t) = 
(T^^(r)^) - {df. In the high-field hmit hx ^ uo, all 
Green functions are strongly peaked at r = 0. U sing the 
adiabatic approximation outlined in Sec. HA, we find 
from Lo[n] [Eq. (|)] (up to Oiuj^/hl)) 



GMr) - 


~ Nh{hx - i<j>) 




(A4) 


G.^r) c: 


.^sgn(r)e^('--^)l^l. 




(A5) 


G,,{r) = 


, 4J 2J(/..-^0)^_f,_ 


-'■0)|t| 


(A6) 



where (j) — (/'(O). Along the classical path, the field cf) 
varies on a timescale I/loq, i.e. slowly on the timescale 
over which the Green functions vanish, which allows us 
to set exp[— (/la; — i(/))|T|] -^ {2/{hx — i(t>))S{T) inG^^ and 
G^^. Because G^^ — ^G^^ the second and third term 
in Eq. ( [A^ ) cancel. To leading order in uj^/hx wc then 
obtain 



Xyyi-r) - (ffMs)^^ (1 - (cos^ ' 



sm' 



0))5(r) 
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^ 2 
(ffA^B)^^e-'^^l^l ^ ig^,)^^Sir), (A7) 



and 



X..(r)c.(g^^f^^^^-'^^l^l^-^-.-^ 



8J 



{9f^Br^S{r). (A8) 



APPENDIX B: TWO-STATE MODEL OF THE 
MODIFIED FW 

As shown in Sees. H and pi for weak tunneling S/h '3> 
1, the low-energy sector of the (modified) FW can be 
described as a two-state model with basis | t) and | |). In 
this approximation, H = -(A/2)[| t)(i | + | i)(T |], where 
the parameters of the original microscopic Hamiltonian 
enter A. 

For any operator O, the transition matrix element be- 
tween \g) and |e) can be evaluated from 



(Bl) 



(e|0|5) = 2((t|0|T)-a|OU) 

+(T|OU>-aiO|T> 



For h^ ^ cDq the state | f) describes a Gaussian prob- 
ability distribution for n in the plane _L B with vari- 
ance (t \cos^<p\ T) = 2J/NhijQ = 1/{S/K) [Sec. |iA|. 
Expanding sin0 ~ ±(1 — cos^ 0/2) for (j) ~ 7r/2 and 
(j) ~ 37r/2, respectively, we obtain 



|(e|sin0|g)| - 2l(t I sin' 
J 



T)-ah 



NhwQ 



(B2) 



The terms (J, | sin(/)| "f) are of order exp[— 5/?i] and hence 
negligible in the weak tunneling regime. Similarly, we 
also find 



|(e|sin3 0|<7)|~l~3 



J 



Nhwn 



(B3) 



APPENDIX C: EFFECTIVE LAGRANGEAN FOR 
AN OPEN AF SPIN CHAIN 

In this appendix we show that an open spin chain with 
an odd number A'^ — 1 of spins s^, i = 2, 3, . . . , A^, can 
be mapped onto the Lagrangcan of a modified FW. For 
simplicity we restrict ourselves to the chain 



Af-l 



N 



h=jj: 



Sj+l 



^h-E^ 



i=2 



N 
i=2 



(CI) 



Again, we assume that the system exhibits AF order, 
and use the staggering s,; = (— l)'sn -f 1 for the spin 
fields, with 1 • n = 0. Spatial variations of the fields n 
and 1 are strongly suppressed in the small system under 
consideration. We hence can proceed as in Sec. C^ and 
carry out the Gaussian integral over 1. We obtain 



-i^chain [hJ 



{N~lf 



[— (h -f- h/jizez — in X n) 



8J(A^-2) 
+ ((h + hAU.e,) ■ nf] - (A - l)k,s 
— s(h • n — i(t>{l -I- cos 9)) 



2 2 



(C2) 



where h^ = —2skz/{N — l)h. Evidently, by appropri- 
ate definition of J and a)oj the Lagrangean of the open 
spin chain can be mapped onto that of the modified FW 
[Eq. (pTl)] with excess spin 5s = s. 
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